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■ ABSTRACT 

| Long-slit spectra have been obtained with the Keck telescope for a sample of 11 

If} , early-type galaxies covering a wide range in luminosity and hence mass. Rotation 

£NJ ■ velocity and velocity dispersions, together with 19 Lick line-strength gradients have 

\ been measured, to, on average, two effective radii. Stellar population models taking 

C****" ■ into account the effect of the non-solar chemical composition have been used to derive 

ages, metallicities and a/Fe abundances along the radius. We find that line-strength 
^ <— | gradients are due, mainly, to variations of the total metallicity with the radius. One 

O galaxy out of 11 shows very strong age gradients, with a young central component, 

while the age gradient for the rest of the sample is very shallow or consistent with 
zero. We also find small variations in the [a/Fe] ratio with radius. Contrary to what 
is expected in simple collapse models, galaxies show both positive and negative [a/Fe] 
| profiles. This rules out a solely inside-out, or outsidc-in, formation mechanism for all 

early-type galaxies. We do not find a correlation between the metallicity and the [a/Fe] 
gradients, and the local metallicity is not correlated with the local velocity dispersion 
k> , for all the galaxies of our sample, which rules out scenarios where the delay in the 

; i ■ onset of the galactic winds is the only mechanism producing the metallicity gradients. 

We found that metallicity gradients are correlated with the shape of the isophotes and 
the central mean age and metallicity of the galaxies, for galaxies younger than ~ 10 
Gyr. We show that the correlation between the gradients and the central values is not 
due to the correlation of the errors and indicates that the same process that shaped 
the gradient, also modified the structural parameters of the galaxies and triggered star 
formation in their centres. This strongly supports the merger scenario for the formation 
of these systems, where the degree of dissipation during those mergers increases as the 
mass of the progenitor galaxies decreases. Finally, we also find a dichotomy in the 
plane grad [a/Fe] -[a/Fe] between galaxies with velocity dispersions below and above 
^200 km s -1 , which requires confirmation with larger samples. 

Key words: Galaxies - galaxies: abundances - galaxies: elliptical and lenticular, cD 
- galaxies: evolution - galaxies: formation - galaxies: kinematic and dynamics 
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The formation and evolution of massive, early- type galax- 
ies constitutes a long-standing and crucial problem in cos- 
mology. In all hierarchical clustering models within a A- 
dominated Cold Dark Matter (ACDM) cosmology, the mas- 
sive early-type galaxies seen now are expected to have 
formed through the merger of smaller galaxies over time (e.g. 
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White & Frenk 1991; Somerville & Primack 1999; Cole et al. 
2000; de Lucia et al. 2006). These models have found strong 
support in the recent studies based on the COMBO- 17 and 
DEEP2 surveys which show that the number density of red 
galaxies has increased since redshift z = 1 (Bell et al. 2004; 
Faber et al. 2005). This is true even for the more massive 
galaxies on the red sequence (although see Cimatti, Daddi 
& Renzini 2006 for another point of view). 

However, most of the studies of the evolution of the Fun- 
damental Plane and colour magnitude relation with redshift 
are compatible with an epoch of formation of the stars in 
these galaxies at z >2 and a passive evolution since then 
(e.g. Kelson et al. 2000; Gebhardt et al. 2003; Fritz et al. 
2005), at least for the brightest galaxies. Fainter galaxies 
may have formed a large percentage of their stars at later 
times (di Serego Alighieri et al. 2005; van der Wei et al. 
2005; Holden et al. 2005). If, as usually assumed in the hi- 
erarchical models, mergers between galaxies trigger star for- 
mation, we would expect elliptical galaxies to show a spread 
in their ages, and, therefore, in their colors, spreading the 
scatter of the color-magnitude and fundamental plane cor- 
relations, contrary to what it is observed. This problem for 
the models may be partially solved if any new stars formed 
are more metal-rich, as both, an increase in age and metal- 
licity produce a reddening of the colors (see., e.g. Faber et al. 
1992; Ferreras et al. 1999). Furthermore, in the last few years 
there have been several studies showing that that mean age 
(as derived with single stellar population models) of more 
massive early-type galaxies is larger than the one in the 
less massive ones (e.g., Caldwell et al. 2003; Nelan et al. 
2005; Sanchez-Blazquez et al. 2006b) which would indicate 
an anti-hierarchical formation of the stars in the galaxies. 
Again, the anti-hierarchical formation of stars can be rec- 
onciled with the hierarchical formation of structures decou- 
pling the epoch of the star formation from the epoch of 
the mass assembly. In this scenario, there is a systematic 
decrease, with increasing mass, of the relative amount of 
dissipation (and, therefore, star formation) experienced by 
the baryonic mass component when they assemble (e.g. de 
Lucia et al. 2006). 

Moreover, a successful scenario of galaxy formation 
must reproduce, not only the photometric properties of the 
galaxies, but also, their structural parameters. Normal and 
low-luminosity elliptical galaxies rotate rapidly, are nearly 
isotropic, show discy-distorted isophotes and cuspy inner 
profiles. On the contrary, giant ellipticals are essentially non- 
rotating, show anisotropic dynamics, core inner profiles, and 
boxy isophotes (Kormendy & Bender 1996). The proper- 
ties of the former can be explained in dissipational mergers, 
while the latter are recovered successfully with dissipation- 
less mergers. 

One way to test this prediction is to study the stellar 
population gradients in early-type galaxies. The stellar pop- 
ulation gradients do not only give information about the 
formation epoch of the stars in the galaxies, but also about 
the formation process, as gas dissipation and mergers af- 
fect the gradients in different ways (e.g. Larson 1975; White 
1980; Mihos & Hernquist 1994; Angeletti & Giannone 2003; 
Kobayashi 2004). The correlation of the derived stellar popu- 
lation parameters (age, metallicity, abundances ratios) with 
other global galaxy properties could offer, in principle, in- 



valuable information to determine how the formation and 
subsequent evolution of early-type galaxies has occurred. 

For example, in the classical models of monolithic col- 
lapse (Eggen et al. 1962; Larson 1974a; Carlberg 1984; Ari- 
moto & Yoshii 1987; Gibson 1997), stars form in all regions 
during the collapse and remain in their orbits with little 
movement inward, whereas gas dissipates into the centre, be- 
ing continuously enriched by the evolving stars. Therefore, 
stars formed in the centre are predicted to be more metal 
rich than those in the outer regions. On the other hand, 
as a consequence of the continuous enrichment by evolving 
stars, the centres are predicted to be less a-enhanced than 
the external parts, unless the collapse is extremely short in 
time. To obtain the hight central values of a/Fe measured in 
early- type galaxies, the duration of the collapse is limited to 
~ 1 Gyr (Arimoto & Yoshii 1987; Thomas 1999), therefore, 
null or very small age gradients are expected within this 
picture. As the the degree of dissipation in these models is 
controlled by the potential well, a strong correlation between 
the metallicity gradient and the mass of the galaxies is also 
predicted. Supernova-driven galactic winds, also may help 
to shape the abundances gradients (Mathews & Baker 1971; 
Larson 1974b; Arimoto & Yoshii 1987; Gibson 1997; Franx 
& Illingworth 1990; Martinelli, Matteucci & Colafrancesco 
1998). These winds, initiated when the energy injected by 
supernovae into the interstellar medium matches that of the 
binding energy of the galaxies, act to evacuate the galaxy of 
gas, thereby eliminating the fuel necessary for star forma- 
tion. As external parts have shallower potential wells, galac- 
tic winds develop earlier than in the central regions, where 
the star formation and the subsequent chemical enrichment 
last longer. This mechanism would also lead to positive a/Fe 
gradients (Martinelli et al. 1998; Pipino, Matteucci & Chi- 
appini 2006). 

Predictions of the resultant metallicity gradient of a 
merger remnant are more complicated, as those depend on 
a large number of free parameters. In general, numerical 
simulations suggest that dissipationless mergers lead to a 
flattening of metallicity gradients (White 1980). However, 
numerical simulations with gas suggests that during the 
merger a significant gas fraction migrates toward the cen- 
tral regions of the merging galaxies resulting in a starburst 
(Barnes & Hernquist 1991). These episodes of star forma- 
tion can also produce metallicity and age gradients (Mihos 
& Hernquist 1994). If star formation is triggered in the cen- 
tre of the galaxy due to the merger, [a/Fe] gradients can be 
both, positive or negative depending on the duration of the 
burst and the original abundance pattern of the gas (see, 
e.g. Thomas, Greggio & Bender 1999; Pipino & Matteucci 
2006). On the other hand, the degree of dissipation and the 
ratio between the masses of the two merging systems pro- 
duce differences in the kinematics, shape of the isophotes 
and other properties of the final remnants (see, e.g. Bekki 
& Shioya 1997; Naab, Jesseit & Burkert 2006). 

Several authors have investigated the variation in the 
colours (e.g. Vader et al. 1988; Peletier et al. 1990; Silva & 
Bothun 1998; Tamura et al. 2000; Saglia et al. 2000; Tamura 
& Ohta 2000,2003; Hinkley & Im 2001; Idiart, Michard & 
de Freitas Pachecho 2002; Menanteau et al. 2004; La Bar- 
bera et al. 2004, 2005; Wu et al. 2005; de Propris et al. 
2005) and line-strength indices (e.g. Spinrad et al. 1971; 
Cohen 1979; Gorgas, Efstathiou & Aragon-Salamanca 1990; 
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Franx & Illingworth 1990; Davidge 1992; Carollo, Danziger, 
& Buson 1993; Davies, Sadler & Peletier 1993; Fisher et al. 
1996; Gonzalez 1993; Mehlert et al. 2003; Proctor et al. 2005; 
Sanchez-Blazquez et al. 2006c; Kuntschner et al. 2006) with 
radius in early- type galaxies, and the results so far appear to 
be contradictory in many cases: For example, Peletier et al. 
(1990) found that, for galaxies brighter than Ms = —21 
the colour gradient gets flatter with luminosity, while the 
contrary is seen for galaxies fainter than this magnitude. 
However, Tamura & Otha (2003) found the opposite trend 
in the relation of their colour gradients with magnitude^. 
Carollo et al. (1993) found a correlation between the Mg2 
gradient and the magnitude, central velocity dispersion and 
mass, but only for faint galaxies, while this correlation has 
not been confirmed by other studies (e.g. Kobayashi & Ari- 
moto 1999). Other uncertain results include the existence of 
a correlation between the line-strength index gradients and 
central values found by some authors (e.g. Gonzalez & Gor- 
gas 1985; Sanchez-Blazquez et al. 2006c, Kuntschner et al. 
2006) but not by others (e.g. Kobayashi & Arimoto 1999; 
Mehlert et al. 2003). 

The lack of agreement between studies is partially due 
to the high signal-to-noise ratio (S/N) necessary to mea- 
sure metallicity gradients with accuracy. In this paper we 
make the first attempt to overcome this limitation, present- 
ing very high quality data allowing for the determination of 
line-strength gradients in a sample of 11 early- type galaxies 
spanning a wide range in luminosity. In Sec. [2] we introduce 
the sample and the observations, together with a description 
of the data reduction. Sec.[3]presents the rotation curves and 
the velocity dispersion (a) profiles. Sec. [4] and [5] describe the 
measurement of the spectroscopic indices and the transfor- 
mation of those into age, metallicity and Q-enhancement 
gradients respectively. In Section |6j we perform a brief anal- 
ysis of the central indices. Sections 17.11 17.21 and 17.31 show 
the analysis of the age, metallicity and [E/Fe] gradients, 
and their correlation with other global parameters of the 
galaxies. Section [8] presents the correlation between the lo- 
cal metallicity and the local a, while in Sec.|9]we discuss our 
results. 



2 THE SAMPLE: OBSERVATIONS AND DATA 
REDUCTION 

The sample consists of 11 galaxies, covering a wide range 
in luminosity, extracted from the field, poor groups, and the 
Virgo cluster. Table[T]summarises the main properties of the 
sample and appendix [X] includes an enhanced discussion of 
each galaxy's properties. 

2.1 Observations and data reduction 

2.2 Observations 

The observations were made using Low Resolution Imaging 
Spectrograph (Oke et al. 1995) in long-slit mode on the Keck 
II telescope. All 11 galaxies were observed on a two-night 
run in 2005 February 8-9. Two 600s exposures were taken 

The colour analysed by Peletier et al. (1990) was U-R while 
Tamura & Otha (2003) present their results using B-R 



for each galaxy except for NGC 4472 where two exposures 
of 600 and 300s were taken. The slit was 175-arcsec long 
and 1.5-arcsec wide, and was oriented along the major axis 
of the galaxies. We use the 400 line mm -1 grism blazed at 
3400 A. This instrumental setup gives a spectral resolution 
of ~ 8.0 A (FWHM) and a dispersion of 1.09 A/pixel. The 
total wavelength coverage is 3110-5617 A. A total of 5 stars 
belonging to both the Lick/IDS stellar library (Gorgas et al. 
1993; Worthey 1994) and MILES library (Sanchez-Blazquez 
et al. 2006d) were observed. These stellar templates though, 
we observed with a different grating than the galaxies and, 
therefore, were only used to measure the seeing during the 
observations and to flux-calibrate the spectra. 

2.3 Data reduction 

The standard data reduction procedures (bias subtraction, 
flat-fielding, cosmic-ray removal, wavelength calibration, sky 
subtraction and flux calibration) were performed with Re- 
duceme (Cardiel 1999). This reduction package allows a par- 
allel treatment of data and error frames and, therefore, pro- 
duces an associate error spectrum for each individual data 
spectrum. After the bias subtraction, the pixel-to-pixel sen- 
sitivity variations were removed (using flat-field exposures 
of a tungsten calibration lamp). 

Prior to the wavelength calibration, arc frames were 
used to correct from geometrical distortions along the spa- 
tial direction (C-distortion) in the images. This correction 
guarantees alignment errors to below 0.1 pixel. Spectra were 
converted to a linear scale using typically 50 arc lines fitted 
by 3th-5th order polynomial with RMS-residuals of 0.1 A. 
The spectra were also corrected for geometrical distortions 
along the spectral direction (S-distortion). To do that we use 
a routine that finds the maximum corresponding to the cen- 
tre of the galaxy as a function of wavelength and fits these 
positions with a low-order polynomial. The spectra are then 
accordingly displaced using a technique that minimizes the 
errors due to the discretization of the signal. The typical 
error in the centring of the object is very small, ~ 0.1 pixel 
(see Cardiel 1999 for details). 

Atmospheric extinction was calculated us- 
ing the extinction curve from the CFHT Bul- 
letin, number 19, pl6 (1998) available through 
http: / / www.jach.hawaii.edu/UKIRT / astronomy/exts.html 
To correct for the effect of interstellar extinction, we used 
the curve of Savage & Mathis (1979), and values of the 
colour excess, E(B — V), from NED. 

Some of the reduction steps, such as sky subtraction, 
can have a significant impact on the derived gradient. Since 
galaxy light levels are usually only a few per cent of the sky 
signal in the outer parts of the galaxies, this process consti- 
tutes one of the most important potential sources of system- 
atic errors (see Cardiel, Gorgas & Aragon-Salamanca 1995). 
Fortunately, our observations were performed in dark time, 
which reduce the variability of the sky during the exposures 
and, therefore, improve the accuracy of the correction. One 
of the main problems to perform a proper sky subtraction 
is that some of the galaxies are large enough to fill a large 
part of the data frames, so that a sky spectrum measured 
at the edges of the frame is contaminated by light from the 
galaxy. In order to avoid this problem, we offset the centre 
of the slit from the galaxy centre. This ensures a lower level 
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Table 1. Sample of galaxies. Type: Morphological classification from the NASA/IPAC Extragalactic database. r cf j: effective radius 
from Burstein et al. (1987), with the exception of NGC 3384, from which the effective radius was extracted from the RC3 catalogue. 
Magnitude: from Hyperleda database (assuming Ho=70km s~ 1 Mpc _1 ). Inner profile types are extracted from Lauer et al. (1995, Ref 1); 
Rest et al. (2001, Ref 2); Ravindranath et al. (2001b, Ref 3); Lauer et al. (2005) (Ref 4); log (v/cr)*: anisotropy parameter calculated 
as described in the text; (a4/a)xl00: the values have been taken from Bender et al. (1989) except for the galaxy NGC 2865, where the 
value is from Reda et al. (2004); Env: environment where the galaxies are located. 



of contamination, at least in the sky measured at one of the 
edges of the frames. The sky level was measured at the end 
of the slit and subtracted from the whole image. To test the 
error in the sky-subtraction we selected the galaxy with a 
larger effective radii of our sample, NGC 4472 and measured 
the number of counts that the galaxy contributed in the sky 
regions. This was done by fitting a de Vaucouleurs profile to 
the surface brightness profile of the galaxy. The contamina- 
tion of the galaxy in the selected region was lower than 1 
per cent. This contamination produce a negligible effect in 
the measured indices. 

Relative flux calibration of the spectra was achieved us- 
ing exposures of standard stars. All the calibration curves 
were averaged and the errors on the indices due to the er- 
rors on the flux calibration were estimated by the differences 
between the indices measured with different curves. 

From each full galaxy frame, a final two-dimensional 
spectrum was created by extracting spectra along the slit, 
binning in the spatial direction the necessary number of ID 
spectra to guarantee a minimum signal-to-noise ratio per 
A (S/N) of 50 in the spectral region of the H/3 index, which 
ensures a maximum relative error for this index of 9 per 
cent (see Cardiel et al. 1998) . Imposing this criteria, we can 
measure line-strength indices, in most galaxies, at distances 
from the centre of r ~2r e //. This S/N ratio at this radius 
has not been achieved previously. 



3 KINEMATIC PROFILES 

Radial velocities and velocity dispersions for each spectrum 
were calculated using MOVEL and OPTEMA algorithms 
described by Gonzalez (1993). The MOVEL algorithm, (im- 
proved version of the classic Fourier quotient method by 
Sargent et al. (1977), is an iterative procedure in which a 
galaxy model is processed in parallel to the galaxy spectrum. 
In this way, a comparison between the input and recovered 
broadening functions for the model allows one to correct the 
galaxy power spectrum from any imperfections of the data 
handling in Fourier space. The main improvement of the 
procedure is introduced through the OPTEMA algorithm, 



which is able to overcome the typical template mismatch 
problem. For each spectrum, 25 star spectra from the MILES 
library (Sanchez-Bl'azquez et al. 2006d) were scaled, shifted 
and broadened according to a first guess of the 7 (mean line 
strength), v (radial velocity) and a (velocity dispersion) pa- 
rameters, having into account the difference in spectral res- 
olution between the stellar and the galaxy spectra. The next 
step is to find the linear combination of the template spec- 
tra that best matches the observed galaxy spectrum. This 
provides a first composite template which is fed into the 
MOVEL algorithm. The output kinematic parameters were 
then used to create an improved composite template and the 
process was iterated until it converged. This iterative ap- 
proach provides an optimal template while simultaneously 
computing the radial velocity and velocity dispersion of the 
galaxy spectrum. In Fig. [1] we show a typical fit between the 
central spectrum of a galaxy and its corresponding optimal 
template corrected with the derived kinematic parameters. 
For each galaxy spectrum, random errors in the derived kine- 
matic parameters were computed by numerical simulations. 
In each simulation, a bootstrapped galaxy spectrum, calcu- 
lated from the error spectra provided by the reduction with 
Reduceme by assuming Gaussian errors, is fed into the algo- 
rithms described above. Errors in the parameters are then 
calculated as the standard deviations of the different solu- 
tions. 

Figure [2] shows the kinematic profiles and velocity dis- 
persion profiles for all the galaxies in the sample. Below we 
analyse the rotation curve and velocity dispersion profile for 
the individual galaxies: 

NGC 1600: Very little rotation and high velocity dis- 
persion along the radius. 

NGC 1700: Clear evidence seen in the profile of the 
presence of a counterotating core, previously found by other 
authors (Statler et al. 1996) and a large rotation at large 
radii. The velocity dispersion profile shows a depression in 
the central parts, which indicates that a more rotationally 
supported component is present at the galaxy centre. 

NGC 2865: This galaxy shows also the presence of a 
kinematically-decoupled core with a radius of ~ 5" and a 
large rotation at larger radii. 
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Figure 1. Top panel: Final fit of the optimal template corrected 
with the kinematic parameters (grey line) to the spectrum of 
the galaxy NGC 4387 (dark line). The bottom panel shows the 
residuals of this fit. The vertical lines indicate the position of 
several typical absorption and emission lines. 



NGC 3377: This galaxy shows large rotation in the cen- 
tre, indicating the possible presence of a central disc. Beyond 
~5 arcsec, the velocity along the radius remains constant. 

NGC 3379: This galaxy shows very little rotation. 

NGC 3384: The rotation curve indicates that the slit 
was not perfectly centred during the observations. The ve- 
locity dispersion profile shows a depression in the centre. 

NGC 4387: The velocity dispersion profile shows a de- 
pression in the centre, which may indicate the presence of a 
stellar disc. 

NGC 4458: The kinematic profile in this galaxy in- 
dicates that the slit was not perfectly centred in the ob- 
ject. Apart from that, the galaxy shows evidence for a 
kinematically-decoupled core. 

NGC 4464: This galaxy shows a rotational curve con- 
sistent with the presence of a gaseous disc (see also Halliday 
et al. 2001) with a scale-length of ~ 4 arcsec. The rotation 
curve starts to decrease beyond this radius. 

NGC 4472: We detect some indication of a 
kinematically-decoupled structure in the centre of the 
galaxy, with very little rotation. 

NGC 4551: This galaxy shows rotation in the central 
parts and then the rotation curve is flat at larger radii. The 
velocity dispersion profile shows a depression in the centre. 

We confirm the presence of kinematically-decoupled 
cores in several galaxies: NGC 1700, NGC 2865, and 
NGC 4458, as previously reported by other studies (Statler 
et al. 1996; Hau et al. 1999; Morelli et al. 2004). Ellipticals 
with kinematically-decoupled cores can originate from dis- 
sipative major galaxy mergers between two spirals (Hern- 
quist & Barnes 1991), a spiral and a elliptical (Franx & 
Illingworth 1988) or by accretion of small satellites via dy- 
namical friction (Kormendy 1984), although the decoupled 
subsystems are often dynamically colder than expected from 
purely stellar satellite accretion (Bender & Surma 1995; 
Franx & Illingworth 1988). Recent studies favour major 
mergers to form kinematically-decoupled structures (Jesseit 



et al. 2006). These kinematically-decoupled structures, how- 
ever, cannot be formed in mergers without dissipation. 

We also found that, for the galaxies NGC 3384 and 
NGC 4458 the slit was slightly offset from the centre of the 
galaxy. In principle, one concern is that this could affect 
the gradients. However, the line-strength indices gradients 
of these two galaxies are compared with the ones derived by 
the SAURON group (Kuntschner et al. 2006) in appendix [C] 
and they are very similar. Therefore, we consider that the 
gradients we are measuring are not very dissimilar to the 
ones we would have obtained with the slit perfectly centred 
in the galaxies. 



4 LINE-STRENGTH INDICES 

In each individual spectrum we measured a total of 15 
Lick/IDS indices (from CNi to Fe5335; Trager et al. 1998), 
the D4000 break (Bruzual 1983) and the 4 higher-order 
Balmer lines H<5a, H<5f, H7.A, H'Jf defined by Worthey & 
Ottaviani (1997). In order to compare the measured indices 
with stellar population models using the fitting functions 
obtained with the Lick/IDS stellar library (Gorgas et al. 
1993; Worthey et al. 1994) the data have to be transformed 
into the Lick/IDS spectrophotometric system. The trans- 
formation into the Lick/IDS system compensates for two 
different effects. The first effect is the instrumental and 
Doppler broadening. Before measuring the indices on the 
spectra, they have to be degraded to the resolution of the 
Lick/IDS library and the indices have to be corrected for 
the internal motions of the stars. The usual method to 
perform this correction is to convolve the stellar spectra 
with different broadening functions of different velocity 
dispersions and use a polynomial expansion to describe the 
variation of a given line index as a function of the input 
velocity dispersion. However, the variation of a given line in- 
dex with velocity dispersion depends on the strength of the 
index itself, (see, for example, appendix B2 in Kuntschner 
2000). Therefore, by assuming a single polynomial we are 
introducing systematic errors that can affect the relation 
of the indices with the velocity dispersion. To avoid these 
problems, we perform the broadening corrections following 
the method described in Kelson et al. (2006) . The corrected 
line strength indices are obtained as follows: 

Icorr = Ig + (It ids — ^S(ct)ot), (1) 

where Ig is the index measured directly on the galaxy 
spectrum, \b(o)oT, the index measured from the optimal 
template spectrum (obtained in the computation of the ve- 
locity dispersions, see Sec[3J) broadened to the instrumental 
resolution and to the velocity dispersion of the galaxy, and 
It ids is the index measured on the optimal template spec- 
trum with a velocity dispersion of zero and convolved to the 
resolution of the Lick/IDS system (Worthey & Ottaviani 
1997). 

The second effect that has to be corrected is the small 
systematic differences that appear when the indices mea- 
sured on the Lick stars and in flux calibrated stars are com- 
pared. These small offsets are partially due to the differ- 
ences in the continuum shape, as Lick/IDS stars are not 
flux-calibrated, although these differences do not explain the 
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Figure 2. Rotation velocity curves and velocity dispersion profiles for our sample of galaxies. 
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H/9 
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0.413 


1.27 


Mgi 
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0.011 


1.27 
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0.011 
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0.012 


2.50 


Mgb 
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0.133 
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2.82 


Fe5270 


-0.112 


0.210 


0.174 


1.47 


Fe5335 


-0.332 


0.201 


0.213 


2.60 



Table 2. Comparison of the Lick/IDS line-strength indices mea- 
sured in this work with Trager et al. (1998). Offset: Mean offsets 
to transform into the Lick system; RMS: root-mean square devia- 
tions of the comparison. RMS(exp): RMS expected by the errors; 
z: ^-parameter to determine the statistical significance of the off- 
set. When z < 1.9 (equivalent to a difference lower than 3cr) we 
consider the offset is not significant and set it as zero. 

offsets of all indices (see., e.g. Worthey & Ottaviani 1997). 
The usual procedure to perform this correction is to observe, 
with the same instrumental configuration as the galaxies, 
stars in common with the Lick/IDS library and compare 
the indices measured in them (after degrading the spectra 
to match the Lick/IDS resolution) with the indices measured 
on the Lick spectra. However, we did not observe stars with 
the same instrumental configuration as the galaxies. There- 
fore, we derived these offsets by comparing the indices of 10 
galaxies in common with Trager et al. (1998), which were 
observed with the same instrumental setup as the Lick/IDS 
library. In order to avoid artificial offsets due to line-strength 
gradients in galaxies, we extracted the spectra in an aper- 
ture of 4" x 1.5", very similar to the one in Trager et al. 
(1998) (4" x 1.4"). 

However, Trager et al. (1998) do not measure the H<5 
and H7 indices. To determine the offsets in the high order 
Balmer lines we compare the indices of 5 galaxies in com- 
mon with Sanchez-Blazquez (2004). We selected the mea- 
surements within an aperture of 2 x 4", which is the most 
similar to the 4 x 1.4" aperture used by Trager et al. A de- 
tailed comparison is shown in Appendix[B] The final offsets 
are listed in Table [2] 

4.1 Emission-line correction 

A large percentage of early-type galaxies show small 
amounts of gas. Some line-strength indices are affected by 
these emission lines. In particular, this is the case of H/3, 
Fe5015 and Mgb. Before measuring line-strength indices this 
emission has to be removed. To measure the flux of the emis- 
sion lines we used the routine GANDALF (Sarzi et al. 2006) . 
This routine fits simultaneously both the stellar spectrum 



and emission lines, treating the emission lines as additional 
Gaussian templates. For the stellar templates, we use syn- 
thetic stellar energy distributions from Vazdekis et al. (2007) 
based on the MILES stellar library (Sanchez-Blazquez et al. 
2006d). Our detection limit is 0.3 A for [OIHU5007 and 
[NIJA5200, and 0.25A for H/3 (see Sarzi et al. 2006 for de- 
tails on how to calculate this threshold). NGC 1700, NGC 
2865, NGC 3377, NGC 3379, NGC 3384 show the presence of 
weak emission lines above our detection limit. The emission 
spectra were subtracted to the galaxy spectra and indices 
measured in the emission-free spectra. 

In order to check the reliability of our line-strength gra- 
dients we have compared the derived line-strength indices 
with the ones derived by other authors, in particular, with 
the studies by Fisher et al. (1996) and Kuntschner et al. 
(2006). The results of this comparison can be found in ap- 
pendix [C] In general, the agreement is very good, although 
none of these studies reaches the same distance from the 
galaxy centres as the present study. 



5 STELLAR POPULATION PARAMETERS 
GRADIENTS 

We derive Single Stellar Population (SSP)-equivalent pa- 
rameters age, [Fe/H] and a-abundance ratios for all the 
galaxies of the sample using the ^-minimization technique 
detailed in Proctor & Sansom (2002) and Proctor et al. 
(2004a,b). The technique involves the comparison of as many 
indices as possible to simple stellar population (SSP) mod- 
els. In the method described by Proctor et al. (2004a,b) 
the indices which significantly deviate by more than 3<r are 
clipped from the fit and x 2 is recalculated. However, in this 
study, we have preferred to obtain the stellar population 
parameters with the same indices for all the galaxies in or- 
der to obtain homogeneous ages and abundances, although 
the results do not change if we clip the deviating indices 
(see below). In particular, we used 19 indices, including all 
Lick/IDS and Worthey & Ottaviani (1997) indices. 

We adopted the SSP models of Thomas, Maraston & 
Bender (2003, TMB03 hereafter) and Thomas, Maraston 

6 Korn (2004, TMK04 hereafter) for the extension to the 
higher order Balmer lines. These models include the effect 
that a variation in the [a/Fe] ratios produce in the Lick in- 
dices. To do that, the authors used the new response func- 
tions by Korn, Maraston & Thomas (2005), as opposed 
to the widely used response functions by Trippico & Bell 
(1995). An important improvement of the Korn et al. (2005) 
response function over the ones by Trippico & Bell (1995) 
is the treatment of the element carbon. Houdashelt et al. 
(2002) pointed out that enhancing the abundance of C by 
+0.3 (which is the factor Tripicco & Bell fitting functions 
use) brings the C/O ratio very close to unity producing 
a carbon star. Worthey (2005) noticed that the response 
functions of Tripicco & Bell (1995) are overestimating the 
sensitivity of C4668 to C by not taking this effect into ac- 
count. The Korn et al. (2005) response functions deal with 
this problem by computing the sensitivity of indices to C 
by increasing the C abundance by only +0.15 dex. The a- 
enhancement is parametrized in these models by [E/Fe]; the 
abundance ratio of the a-elements O, Ne, Mg, Si, S, Ar, 
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Error bars represent the scatter in the deviation. 



Ca and Ti plus the elements N and Na. (see the original 
references for details). 

Residuals to the best-fitting (observed value minus best- 
fitting value expressed in terms of the index errors) for all 
the galaxies are summarised in Fig. [3] The points represent 
the average deviation from the best-fitting values while the 
errors bars show the RMS for all the points along the radius. 
As can be seen, the mean differences between the fitted and 
the measured errors are very small, confirming the reliability 
of the fitting procedure. 

Figures [4][5]and[6] show the age, metallicity and [E/Fe] 
profiles for all the galaxies in our sample. We have to remind 
the reader that those values of age, metallicity and [E/Fe] 
are SSP-equivalent parameters. If galaxies experience a star 
formation history more complicated than a single, instanta- 
neous burst, those values have to be carefully interpreted: 
for example, if a galaxy has experience two burst of star for- 
mation separated in time, the SSP-age would be biased to- 
ward the age of the youngest stars, while the SSP-metallicity 
would be closer to the metallicity of the old population (al- 
though this depends, of course, of the mass fraction of the 
bursts, see, e.g. Serra & Trager 2006 for a detailed analysis). 

To quantify the profiles we perform a linear fit weighting 
with the errors in the y-direction and excluding the points 
from the centre within the diameter of the seeing disc to 
avoid seeing effects. The extent of the seeing disc (FWHM 
of the point spread function) was taken from the measure- 
ments performed by the Mauna Kea seeing monitor, but 
the measurements were checked by measuring the FWHM 
of the Point Spread Functions on observed standard stars. 
The points excluded from the fit are represented with open 
symbols. For completeness, we have performed also a linear- 
fit on the age profile of the galaxy NGC 3377, despite this 
profile being clearly non-linear. The errors in the gradients 
were calculated as the unbiased standard deviation from the 



fit. Table [3] shows the derived gradients of age, [Z/H] and 
[E/Fe] . 



5.1 Checking the reliability of the stellar 
population gradients 

One of the major problems of using a large number of in- 
dices to derive stellar population parameters by \ 2 mini- 
mization is related to the zero-point of the models. Although 
we have added suitable offsets to transform our data into the 
Lick/IDS spectrophotometric system, the models have prob- 
lems reproducing some of the measured indices for Galactic 
globular clusters (e.g. TMB03). The same problem has been 
observed in elliptical galaxies where it has been also shown 
that the age and metallicity depend on the index used to 
measure those parameters (e.g. Peletier et al. 1990; Faber 
et al. 1992; Vazdekis et al. 2001; TMB03, Sanchez-Blazquez 
et al. 2006b) 

Most of these differences, but not all, can be attributed 
to peculiar abundances ratios in the Galactic globular clus- 
ters and elliptical galaxies (TMB03). In fact, TMB03 showed 
that part of the problem is solved with the inclusion of non- 
solar [a/Fe] abundances in the models but, certainly, this is 
still an issue when one wishes to use a large number of in- 
dices with different sensitivities to different chemical species 
as (1) some other elements apart from the ones considered 
in the models may have abundances ratios different from so- 
lar, and (2) because it may be the case that not all the a 
elements are equally enhanced, as seen, for example, in the 
Galactic bulge (Lecureur et al. 2007). 

Kelson et al. (2006) treated this problem by recalibrat- 
ing the stellar population models. To do that, they used 
a reference galaxy and assigned to it an approximate age, 
metallicity and [a/Fe]. They derived offsets in all the indices 
so the final indices coincide with the values given by the 
models for the selected stellar population parameters. They 
applied these offsets to all the galaxies in their sample. This 
study was the first to take the effect of the zero-point into ac- 
count in the derivation of ages and metallicities using a large 
number of indices. However, their approach assumes that the 
offsets are the same for all the galaxies, while, if the offset is 
due to differences, for example, in the chemical abundance 
ratios not considered in the models, this assumption presum- 
ably would not be true, as the chemical composition is dif- 
ferent for galaxies of different masses. Furthermore, if galaxy 
stellar populations are not suitably described by means of 
SSP models, the age and metallicities measured with differ- 
ent indicators might be different, simply because the weight 
of different stars to different regions of the spectra is different 
(e.g. Sanchez-Blazquez et al. 2006b; Schiavon et al. 2006). 
The only way to deal with this is to use a deconstruction 
method to analyse the whole star formation history of the 
galaxy (e.g., Panter, Heavens & Jimenez 2003; Ocvirk et al. 
2006a,b). This will be done in a future paper. Instead, in 
the present work, we check the robustness of our results by 
comparing our gradients with the ones obtained using differ- 
ent indices and models. In particular, we will compare the 
age, metallicity and [E/Fe] gradients obtained as described 
above with the ones obtained as follow: (1) Using a set of 
three indices: Fe4383, H/3 and Mgb; (2) Using 19 indices, but 
clipping the indices that deviate more than 3-er from the fit; 
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Figure 4. SSP-equivalent-age profiles for the sample of galaxies calculated with TMB03 models using a set of 19 different Lick-indices. 
Dashed lines indicate the extent of the seeing disc. Solid lines represent a linear-fit to the data weighting with errors in the y-direction. 
Open symbols have been excluded from the fit (see text for details). 



(3) Using the models by Vazdekis et al. (20070 and a the 
indices Fe4383, H/3 and Mgb. The last option is being shown 
as an extreme case of a completely different technique and 
model and serves a purpose in showing the robustness of the 
results against the technique and models employed. 

(1) We selected H/3, Fe4383 and Mgb because the mod- 
els do an excellent job fitting these indices to globular clus- 
ters measurements and also because they have the highest 
sensitivities to the parameters we want to measure, i.e., age, 
Fe, and Mg abundances respectively. We perform a linear fit 
to the final values with radius in the same way as described 
above. The gradients are shown in Table [3] 

(2) Although we decided to use the same indices for all 
the galaxies and along the radius to avoid any possible bias 
in the result due to the different sensitivity of the different 
indices to different chemical compositions, here we compare 



2 availables at http:/ /www. ucm.es/info/Astrof/users/pat/models. html 



with the results obtained when a 3-a clipping is applied dur- 
ing the x 2 minimization as in the original method described 
by Proctor & Sanson (2002). 

(3) Finally, we also calculate the stellar population pa- 
rameters along the radius using the new stellar population 
models by Vazdekis et al. (2007). These models are built for 
solar chemical composition and they are only calibrated for 
variations of age and [Z/H]. To derive [E/Fe] values, we fol- 
lowed the same procedure as in Trager et al. (2000a, TOOa 
hereafter), with the only difference that we use the new fit- 
ting functions of Korn et al. (2005) instead of the Trippico 
& Bell (1995) ones used by TOOa. We include in the en- 
hanced group ([X/H]=+0.3) the elements C, N, O, Mg, Fe, 
Ca, Na, Si, Cr, and Ti, while the elements Cr, Mn, Co, Ni, 
Cu and Zn are depressed by [X/H] =— 0.3. We interpolate the 
model grid at intervals of 0.1 Gyr in age and 0.05 in metal- 
licity. Then we applied the response functions by Korn et al. 
(2005) to obtain an interval step of 0.05 in [E/Fe]. Stellar 
population equivalent parameters were derived by choosing 
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Figure 5. SSP-equivalent-[Z/H] profiles for the galaxies of the sample calculated with TMB03 models using a set of 19 different line- 
strength indices (see text for details). Solid lines represent a linear-fit to the data weighting with the errors in the y-direction. Only solid 
points have been fitted. Dashed lines indicate the extent of the seeing disc. 



the best-fitting age, [Z/H] and [E/Fe] for the indices H/3, 
Fe4383 and Mgb. 

Fig. [7] represents the comparison of the gradients calcu- 
lated with all the 19 indices as described at the beginning 
of the section with the methods (1), (2) and (3) described 
above. Table|4]shows the mean offset and residual dispersion 
of the comparison, calculated as the quadratic difference be- 
tween the measured dispersion and the one expected from 
the errors. 

As can be seen, significant offsets do not exist between 
the gradients calculated with a different number of indices. 
We do not even find significant offsets between the gradi- 
ents calculated with different models and different number 
of indices. The RMS dispersion in the comparison of the 
methods with and without clipping is compatible with the 
errors, however, the dispersion among the 1:1 relation in 
the comparison of the methods (1) and (3) with the method 
employing 19 indices is larger than the one expected by the 
errors. This should prevent us drawing conclusions that are 



not confirmed using the gradients measured with all the dif- 
ferent methods. 

In the rest of the paper we use the gradients derived 
with 19 indices without clipping because the errors on the 
parameters are smaller. We have checked, though, that none 
of the conclusions would change if we used the ones calcu- 
lated with any of the other methods described above. 



6 CENTRAL VALUES 

In order to compare the gradients with the central values, we 
also extracted, for each galaxy, a spectrum within an aper- 
ture of 1.5" x r e f/8 and used the corresponding line-strength 
indices to derive stellar population parameters. Fig.|5]shows 
the classical relations between the Mgb and the central ve- 
locity dispersion for our sample of galaxies. Table [DTI in ap- 
pendix , shows all the central indices and the central velocity 
dispersion. The values of the mean SSP-age, -metallicity and 
-[E/Fe], as derived with 19 indices are reported in Table [5] 
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Figure 6. SSP-equivalent-[E/Fe] profiles for the galaxies in the sample calculated with TMB03 models using a set of 19 different indices. 
The solid lines represent linear-fits to the solid points weighting with the errors in the ^/-direction. The open symbols are those excluded 
from the fit. Dashed lines indicate the extent of the seeing disc. 



Fig. [9] shows the relation of the age, metallicity and 
[E/Fe] with the central velocity dispersion. We have sepa- 
rated the galaxies -with different symbols- as a function of 
their central ages. As can be seen, contrary to what happens 
when the indices are compared with this parameter, the re- 
lations appear to have considerable dispersion. We analyse 
separately the three relations: 

[Z/H] vs. a: The relations between the Mg2 index 
and a and the color-magnitude relations have been inter- 
preted, classically, as a mass-metallicity relation for early- 
type galaxies (e.g. Kodama & Arimoto 1997), in the sense 
that more massive galaxies are also more metal rich. How- 
ever, it is clear from Fig. [9] that we do not detect a corre- 
lation between the metallicity and the velocity dispersion. 
Old galaxies (with ages > 8 Gyr) define a mass-metallicity 
sequence, but galaxies with intermediate- or young- ages, lie 
above the relation in such a way that the deviation from the 



relation correlates with the light-weighted mean age. The 
same result was obtained by j0rgensen (1999) and Trager 
et al. (2000b). Trager et al. (2000b) showed that galaxies 
are in a plane within the space of age, [Z/H] and a, where, 
at a given a, galaxies more metal rich are also younger. In 
the figure we have plotted the best-linear-relation found by 
Nelan et al. (2005) and the RMS dispersion found by those 
authors. As can be seen, the oldest galaxies of our sample 
(with ages>8 Gyr) follow a very similar relation than the 
one found by these authors, but the younger galaxies devi- 
ate systematically from the relation. As several authors have 
noted before, the fact that, at a given a, younger galaxies 
are also more metal rich has very important implications 
as, if real, the tightness of the colour-magnitude relation 
should not be interpreted as a low dispersion in the mean 
stellar ages of galaxies. This is because the effect that a re- 
juvenation of the stellar population has in the line-strength 
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189 ±0.080 


-0 


.021 ±0.051 





.020 ± 0.045 



Table 3. Age, [Z/H] and [E/Fe] gradients for the 11 galaxies in our sample calculated with a x 2 minimization with different number of 
indices and different models. See text for details. 



TMB03 models TMB03 models V07 models 

3 indices 19 indices (clipping) 3 indices 

offset RMS RMS RMS offset 



offset 


RMS 


RMS 


RMS 






(exp) 


(res) 


0.002 


0.03 


0.02 


0.022 


-0.043 


0.07 


0.05 


0.048 


-0.009 


0.05 


0.04 


0.030 



grad age 0.002 0.03 0.02 0.022 0.000 0.000 0.000 0.000 0.007 

grad [Z/H] -0.043 0.07 0.05 0.048 -0.006 0.049 0.051 0.000 -0.014 

0.014 0.030 0.035 0.000 -0.053 



RMS 


RMS 


RMS 




(exp) 


(res) 


0.000 


0.000 


0.000 


0.049 


0.051 


0.000 


0.030 


0.035 


0.000 



RMS 


RMS 


RMS 




(exp) 


(res) 


0.197 


0.122 


0.155 


0.080 


0.077 


0.022 


0.059 


0.052 


0.028 



Table 4. Comparison of the SSP-equivalent-parameter gradients derived using different methods: Mean offset, RMS-dispersion, RMS- 
dispersion expected by errors and RMS-residual (not explained by the errors) in the comparison of the gradients measured with a \ 2 
minimization and 19 indices without clipping and the following methods: (1) TMB03 (3 indices): Using TMB03 models and only the 
indices Fe4383, H/3 and Mgb; (2) TMB03 (19 indices, clipping): Using TMB03 models and a \ 2 minimization technique, clipping the 
indices that deviate more than 3— cr from the final fit; (3) V07 models (3 indices): Using Vazdekis et al. (2007) models and only the 
indices H/3, Fe4383 and Mgb indices. 



indices and colour is very similar to a decrease in metallicity 
- the so-called age-metallicity degeneracy (see, e.g. Worthey 
1994). This should be taken into account, also, when inter- 
preting the evolution of the colour-magnitude diagram with 
redshift. 

The result shown here and in the studies mentioned 
above challenges the interpretation of the color-magnitude 
relation as a purely mass-metallicity relation, as usually as- 
sumed in the literature, as galaxies with the same mass can 
show a large spread in their metallicities. 

[E/Fe] vs a: Several authors have found the existence of 
a positive correlation between the degree of a-enhancement 
in the central parts of early-type galaxies and a (Faber et al. 
1992; Worthey, Faber & Gonzalez 1992; Kuntschner 2000, 



J0rgensen 1999; Trager et al. 2000b; Terlevich & Forbes 
2002; Thomas et al. 2005; Nelan et al. 2005). a-elements and 
Fe-peak elements are released into the interstellar medium 
by different stars and on different time-scales. Therefore, it 
is common in the literature to use [E/Fe] as a cosmic clock 
to quantify the duration of the star formation. The reported 
relation between the ratio [E/Fe] and a would imply, in this 
context, that more massive galaxies formed their stars in 
shorter time-scales. The relation between the [E/Fe] and the 
central velocity dispersion is, however, not very clear in our 
data (a non-parametric Spearman rank order test gives a 
probability of no correlation of ~5 per cent, which is less 
than 2-ct significance). The most massive galaxies are the 
ones with higher [E/Fe], but there is considerable scatter 
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Figure 7. Comparison between the gradients of SSP-equivalent parameters calculated with the TMB03 models using all 19 indices and 
the ones obtained with the following methods: (1) Using TMB03 models and only the indices Fe4383, H/3 and Mgb (upper panels); (2) 
Using TMB03 models and 111 9 indices, but clipping the indices deviating more than 3-er from the best-fitting (mid panels); (3) Using 
Vazdekis et al. (2007) models (combined with the TOO method to derive [E/Fe]) and the indices Fe4383, H/3 and Mgb (lower panels). 



at lower a. We have also plotted, in the figure, the relation 
found by Nelan et al. (2005) using ~500 galaxies. As can be 
seen, our galaxies follow the same relation, and the lack of 
correlation may be due to the small sample size. 



Finally, we present in the bottom panel of Fig. the 
relation between the central age and the central velocity 
dispersion. Contrary to what has been reported in many 
studies (e.g. Caldwell et al. 2003; Nelan et al. 2005; Bernardi 
et al. 2006; Sanchez-Blazquez et al. 2006b) we do not find 
any relation between these two parameters. Although we do 
not find young galaxies with a > 200 km s , we find a very 
large dispersion in the ages of galaxies with a < 200 km s _1 . 
In the figure we also show the relation found by Nelan et al. 
(2005). It can be seen that the youngest galaxies in our sam- 
ple clearly deviate from the relation found by those authors, 
although the rest of the galaxies are consistent with it. 



7 GRADIENTS OF AGE, [FE/H] AND [E/FE] 

The variation of the physical properties of galaxies with ra- 
dius prove invaluable information for constraining the pro- 
cesses of galaxy formation and evolution. In the following 
sections we will analyse the relation between the stellar pop- 
ulation parameters and other properties of the galaxies try- 
ing to understand the mechanism dominating the formation 
of the measured gradients. 

7.1 Age gradients 

We find that for 10 of the 11 galaxies, age gradients are 
compatible with zero-slope. NGC 1700 show a significant age 
gradient when this is measured using 19 indices, but this is 
not true if only 3 indices are employed, independently of the 
model used. NGC 3377 seems to show an age gradient, but 
the age profile in this galaxy does not behave as a power- 
law, and, therefore, it does not make much sense to perform 
a linear-fit. Several authors have analysed age gradients in 
samples of early-type galaxies but the results are still not 



14 P. Sanchez-Bldzquez et al. 



Galaxy Age (Gyr) [Z/H] [E/Fe] 



NGC 1600 


13.5 ±0.1 


0.326 


± 


0.003 


0.357 


± 


0.004 


NGC 1700 


2.9 ± 0.03 


0.627 


± 


0.006 


0.186 


± 


0.002 


NGC 2865 


1.0 ±0.01 


0.523 


± 


0.005 


0.133 


± 


0.001 


NGC 3377 


5.5 ± 0.05 


0.343 


± 


0.003 


0.277 


± 


0.003 


NGC 3379 


14.2 ±0.1 


0.244 


± 


0.002 


0.288 


± 


0.003 


NGC 3384 


3.06 ± 0.03 


0.591 


± 


0.006 


0.186 


± 


0.002 


NGC 4387 


8.88 ± 0.09 


0.123 


± 


0.001 


0.154 


± 


0.001 


NGC 4458 


11.7 ±0.1 


0.021 


± 


0.001 


0.241 


± 


0.002 


NGC 4464 


13.5 ±0.1 


0.067 


± 


0.001 


0.256 


± 


0.003 


NGC 4472 


11.3 ±0.1 


0.352 


± 


0.003 


0.285 


± 


0.003 


NGC 4551 


6.47 ±0.06 


0.247 


± 


0.002 


0.183 


± 


0.002 



Table 5. Central values (within an aperture of 1.5" X r c f/8) for the Age, [Z/H] and [E/Fc] measured in our sample of galaxies. 
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Figure 8. Relation between the Mgb index, measured in a aper- 
ture of 1.5" X r e ff/8 versus the central velocity dispersion inside 
the same aperture. 



conclusive. For example, Hinkley & Im 2001; Mehlert et al. 
(2003) and Wu et al. (2005) did not find any overall age 
gradient in their samples. Some other authors (e.g. Silva & 
Elston 1994; Gonzalez 1993; Tantalo, Chiosi & Bressan 1998; 
Sanchcz-Blazquez et al. 2006c), however, have detected non- 
zero age gradients in their respective samples. The typical 
values found by those authors is ~ 20 per cent variation, 
with the central parts younger than the external parts. 



7.2 Metallicity gradients 

The mean [Z/H] gradient of our sample is A [Z/H] /logr — 
—0.306 with an RMS dispersion of 0.133. This value sug- 
gests a mean reduction of metallicity in elliptical galaxies of 
more than 50 per cent per decade of variation in radius. This 
value is compatible with the values derived by other stud- 
ies, e.g. A[Z/H]/logr=-0.23 ± 0.09 (Gorgas et al. 1990); 
A[Z/H]/logr=-0.25 ± 0.1 (Fisher et al. 1995). In princi- 
ple, the strength of the metallicity gradient is related to 
the merging history of the galaxies, as, while dissipational 
processes tend to create stronger gradients, mergers be- 
tween galaxies destroy these gradients (e.g., Mihos & Hern- 
quist 1994; Kobayashi 2004) The mean metallicity gradients 
for non-merger and merger galaxies derived by Kobayashi 
(2004) are A[Z/H]/Alogr ~ -0.30 and -0.22, respectively. 
As can be seen, galaxies in our sample are compatible with 
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Figure 9. Relation between the central metallicity, a- 
enhancement and age with the central velocity dispersion. Dif- 
ferent symbols in the two panels indicate different central ages as 
indicated in the inset. Solid and dashed lines indicate the best- 
linear-fitting obtained by Nelan et al. (2005) and the 1— <r errors 
obtained by those authors respectively. 
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both of these values. The direct comparison between the 
metallicity gradient obtained with single stellar population 
models and numerical simulations is, however, difficult, due 
to the manner by which the results from numerical simu- 
lations are transformed to the observational plane. Another 
way to determine the evolutionary paths of early-type galax- 
ies is to study the relation between the metallicity gradients 
and other global properties of these systems, as different 
physical processes are expected to lead to different correla- 
tions. For example, dissipational processes are believed to 
create steeper gradients in more massive galaxies (Carlberg 
1984; Bekki & Shioya 1999), although this is sensitive to 
the adopted feedback prescription in the simulations (e.g. 
Bekki & Shioya 1998; Bekki & Shioya 1999). Dissipation- 
less mergers of galaxies, on the contrary, are expected to 
produce some dilution of the gradients in galaxies (White 
1980) deleting or producing an inverse correlation among 
stellar population gradients and luminosity. 

Fig. [10] shows the correlation of the metallicity gradi- 
ent with the velocity dispersion for our sample of galax- 
ies. Although the sample is not very large, we confirm the 
lack of correlation previously noted by other authors using 
line-strength indices (e.g. Gorgas et al. 1990; Davidge 1991, 
1992; Davies et al. 1993, Mehlert et al. 2003). Galaxies with 
steeper gradients are not the most massive of our sample 
but the ones with a ~ 200 km s _1 . 

A similar trend was found by Kormendy & Djorgovski 
(1989) between the color gradients and Ms in a sample ob- 
tained by combining data with different quality from Vader 
et al. (1988), Franx et al. (1989) and Peletier et al. (1989). 
Carollo et al. (1993) compared the Lick index Mg2 gradi- 
ent and both, the mass, and the luminosity (derived from 
the Fundamental plane), and they concluded that there is 
a change in the slope of the trends approximately at the 
same mass as found here. They also found a significant cor- 
relation between the gradients and mass for galaxies with 
Mtot ^ 10 1 Mq, where M to t represents the total mass of the 
galaxy. Our sample is too small to perform reliable statisti- 
cal tests to study the significance of the correlations in the 
two magnitude ranges. Confirming the presence of a slope 
change in the grad[Z/H]— a plane at Ms ~ —21.5 will re- 
quire a larger sample of comparable quality to that presented 
here. However, there is a significant degree of evidence now 
which indicates that the metallicity gradient is not constant 
with the mass of the galaxies, but gets steeper for galaxies 
around this magnitude. 

As we said earlier, the correlation between the metal- 
licity gradient and the mass of the galaxies in the context 
of mergers of galaxies depends on the degree of dissipation 
in the merging. A positive correlation, as suggested by some 
authors (e.g. Carollo et al. 1993; this study) for galaxies with 
Ms > —20.5, is expected in mergers with gas (Bekki & Sh- 
ioya 1999), while the opposite trend (as suggested by some 
studies, e.g., Vader et al. 1998; Franx 1988) is expected in 
dissipationless mergers, assuming that more massive galax- 
ies have suffered more mergers (as predicted by hierarchical 
models of galaxy formation) . Therefore, if the trends are con- 
firmed, they could be explained assuming a decrease, with 
the mass of the galaxy, of the degree of dissipation during 
the last major merger event. 

Several studies (Carollo et al. 1993; Gonzalez & Gor- 
gas 1995) have found that galaxies with stronger central 
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Figure 10. Relation between the metallicity gradient and the 
central velocity dispersion. Different symbols indicate different 
shapes of the inner profile as indicated in the inset. 



Mg2 indices showed, also, steeper Mg2 gradients. Although 
other authors (Kobayashi & Arimoto 1999; Mehlert et al. 
2003) have failed to find this correlation, it has been re- 
cently confirmed by some studies (Sanchez-Blazquez et al. 
2006c; Kuntschner et al. 2006). However, Kuntschner et al. 
(2006) claim that the correlation is driven by SO galaxies or, 
in particular, by young galaxies. If the central metallicity 
values are correlated with the metallicity gradient, it would 
imply that the global metallicities of ellipticals were more 
similar than the central ones. Unless such aperture effects 
are taken into account, the impact upon the interpretation 
of scaling relations - such as the redshift evolution of the 
colour-magnitude relation - could be significant. 

In Fig. |TT| we have plotted this relation, but using the 
metallicity instead of the Mg2 index. We have separated, 
with different symbols, galaxies in different ranges of cen- 
tral ages. There exists a trend for which galaxies with higher 
central metallicities also show steeper gradients in this pa- 
rameter, although there is considerable scatter among the 
relation. A non-parametric Spearman rank-order test gives 
a level of significance lower than 0.05 (or a probability of cor- 
relation of 95 per cent) , which is the limit to claim a correla- 
tion. In principle, this correlation could be the consequence 
of the correlation of the errors, as an increase in metallicity 
in the centres would produce also an increase in the metal- 
licity gradient. To check this possibility we have represented, 
in Fig.[T2]the relation between our metallicity gradients and 
the central ages from other authors. In particular, we have 
used the values obtained by Terlevich & Forbes (2000) using 
Worthey (1994) stellar population models and the combina- 
tion of indices Mgb, <Fe> and H/3 collected by different 
studies. As can be seen, the trend is still present even us- 
ing completely independent values, which argues against the 
correlation of the errors being the only reason for its exis- 
tence. 

To explore the relation between the metallicity gradi- 
ent and the central age of the galaxies we have represented, 
in Fig. 1131 these two magnitudes. And, as can be seen, the 
correlation is obvious for galaxies with central ages less than 
10 Gyr. For the whole sample, a rank-order Spearman cor- 
relation test gives now a probability of non-correlation lower 
than 1 per cent. We have also separated the galaxies with 
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Figure 11. Relation between the metallicity gradient and the 
central metallicity. Different symbols show galaxies younger than 
5 Gyr (filled circles), between 5 and 8 Gyr (stars) and older than 
8 Gyr (squares). 
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Figure 13. Relation between the metallicity gradient and the 
central age. Different symbols represent galaxies with metallicity 
higher (triangles) and lower (squares) than [Z/H]=+0.2 as indi- 
cated in the inset. 
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Figure 12. Correlation between the metallicity gradient mea- 
sured in this work and the independent central metallicity values 
extracted from Terlevich & Duncan (2000). 




central age (Gyr) (TF00) 

Figure 14. Relation between the metallicity gradient calculated 
in this study and the independent central mean age of the galaxies 
extracted from Terlevich &; Forbes (2000). 



central metallicities higher (stars) and lower (circles) than 
[Z/H]— +0.2. In general, galaxies showing steeper metallic- 
ity gradients are more metal rich and younger while galaxies 
with shallower metallicity gradients are older. For galaxies 
older than 10 Gyr, however, half of the sample have low 
metallicities in their centre, while half of the sample show 
metallicities higher than [Z/H]=+0.2. (see Table EJ. The 
galaxies with larger metallicity are the most massive galaxies 
of our sample, with a > 200 km s _1 . The correlation between 
the metallicity gradient and the central age has been pre- 
viously found by Sanchez-Blazquez et al. (2006c) (although 
the methodology is very different to this work and the cor- 
relation is only found when the metallicity is measured with 
some indicators). As we did before, we checked the possibil- 
ity that the correlation was artificial, due to the correlation 
of the errors. To do that we plotted the metallicity gradients 
against the completely independent ages extracted from the 
Terlevich & Forbes (2000) catalogue. This can be seen in 
Fig. Q3] It is clear, as before, that the correlation is not an 
exclusive consequence of the correlation of the errors. 

Finally, we study the relation between the metallicity 



gradient and the rotation of the galaxies and the shape of 
their isophotes. The decreasing rotational support and tran- 
sition from discy to boxy isophotal shapes with increasing 
stellar mass in elliptical galaxies suggests an increasing frac- 
tion of dissipationless mergers in the growth of the most 
massive elliptical galaxies (e.g. Bender et al. 1992; Kor- 
mendy & Bender 1996; Faber et al. 1997; Naab, Khochfar & 
Burkert 2006). Fig. ll5l shows the relation between the metal- 
licity gradient and the anisotropy parameter (log(w/<r)*) as 
defined in Bender (1990) normalized by the value expected 
for an isotropic oblate galaxy flattened by rotation. To derive 
this parameter we followed the procedure described in Pe- 
draz et al. (2002). A conservative estimate of the maximum 
rotation velocity (v) was computed as the error-weighted 
mean of the two data pairs with the highest rotational veloc- 
ities. To compute the mean velocity dispersion we co-added 
all the individual spectra with radii between the seeing limit 
and the effective radius. Prior to this, we shifted all the spec- 
tra to rest-frame using the rotation curves. The final values 
for the sample of galaxies presented here are listed in Ta- 
ble [l] A non-parametric Spearman rank-order test gives a 
significant correlation with a level of significance lower than 
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Figure 15. Relation between the metallicity gradient and the 
anisotropy parameter normalized to an oblate rotator (Bender 
1990). 
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Figure 16. Relation between the metallicity gradient and the 
parameter (a4/a)x 100, extracted from Bender et al. (1989). 



0.05. However, we are aware that this correlation is driven 
mainly by two galaxies (NGC 1600 and NGC 2865) and that 
larger samples would be necessary to confirm its existence. 

Fig. [16] shows the relation between the metallicity gra- 
dient and the parameter (a4/a)x 100, which measured the 
deviation of the shape of the isophotes from a perfect el- 
lipse. Positive values of a4 indicate discy isophotes while 
a4<0 indicates that the isophotes have a boxy shape. Most 
of the values were taken from Bender et al. (1989), who de- 
fined this parameter as the peak value of a4/a in the case 
of peaked profiles and the value at 1 effective radius in the 
case of monotonically increasing or decreasing profiles. For 
the galaxy NGC 2865, the value was extracted from Reda 
et al. (2004). We could not find references (measured in the 
same way) for three of our galaxies: NGC 3384, NGC 4458 
and NGC 4464 and, therefore, they are not included in the 
Figure. For the rest of the galaxies, there exists a strong cor- 
relation between grad[Z/H] and (a4/a)xl00. Discy galaxies 
show stronger gradients and the strength of the gradient gets 
lower with the boxiness of the isophotes. This strong correla- 
tion is surprising as the a4 parameter measured in simulated 
galaxies depends on the projection effects, and therefore, the 



same galaxy can have boxy and discy isophotes depending 
of the viewing angle (Stiavelli, Londrillo & Messina 1991; 
Governato, Reduzzi & Rampazzo 1993; Governato et al. 
1993; Heyl, Hernquist & Spergel 1994; Lima-Neto & Combes 
1995; Gibson et al. 2006). Bekki & Shioya (1997) performed 
numerical simulations of mergers between gas-rich galaxies, 
studying the effect of star formation on the structural pa- 
rameters of the remnant. They found that the rapidity of gas 
consumption by star formation greatly affects the isophotal 
shape of the merger remnant. Mergers with gradual star for- 
mation are more likely to form elliptical galaxies with discy 
isophotes, while those where the star formation is more rapid 
are more likely to form boxy ellipticals (although this de- 
pends on the viewing angle and, therefore, these scenario 
can lead to galaxies that can be seen as discy too). This 
scenario could explain the relation between the metallicity 
gradient and the shape of the isophotes found here. Fur- 
thermore, if the metallicity gradients are correlated with the 
shape of the isophotes and also with the central values of age 
and metallicity, this implies that the young ages observed in 
many early-type galaxies are not the consequence of a frost- 
ing population that form, for example, due to the accretion 
of small galaxies, but that the physical process that pro- 
duced the recent star formation is related with the process 
shaping the gradients in the galaxies. This could happen, for 
example, in a scenario of mergers between galaxies. 



7.3 Chemical abundances gradients 

The chemical abundance ratios variations with radius give 
information about the time-scales of the star formation 
within the galaxiefl Fig. [6] shows the [E/Fe] gradients for 
our sample of galaxies. As can be seen in the figure, galaxies 
show a broad variety in the slope of the [E/Fe] profiles. In 
particular, contrary to what it would be expected by mod- 
els of dissipative collapse, some of the galaxies show negative 
[E/Fe] gradients, indicating, if interpreted assuming a sin- 
gle episode of star formation, more extended star formation 
histories in the external parts. 

Several theoretical works have studied the [Mg/Fe] gra- 
dients using chemical and cosmological chemodynamical 
evolution cores. Martinelli et al. (1998) and Pipino, Mat- 
teucci & Chiappini (2006) modeled the scenario suggested 
by Franx & Illingworth (1990) where metallicity gradients 
are the consequence of the time delay in the development 
of galactic winds between the central and external parts of 
the galaxies. In this scenario stars in the outermost regions 
form earlier and faster than the ones in the centre and, there- 
fore, a natural outcome from this model is an increase of the 
[Mg/Fe] ratio with radius. In particular, the recent work by 
Pipino et al. (2006) predicts a value of grad[E/Fe]~+0.2- 
0.3 dex, which reproduces the observations by Mendez et al. 
(2005) but it is steeper than any of the gradients measured 
here. In this scenario, the duration of the star formation is 
the only parameter controlling the local metallicity. There- 
fore, we would expect to find a correlation between the [Z/H] 
and the [E/Fe] gradients (under the assumption that [E/Fe] 



4 This is true if we assume that there is not variations of the IMF 
along the radius and that the feedback processes blow away all 
the elements with the same efficiency. 
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Figure 17. Relation between the [E/Fe] gradient and the metal- 
licity gradient for our sample of galaxies. 




2 2.2 2.4 

log a (km/s) 

Figure 18. Relation between the [E/Fe] gradient and the central 
velocity dispersion of the galaxies. 



is a good measure of the time-scales of the star formation) . 
Fig. [TTjshows the relation between both gradients. As can 
be seen, we do not detect any significant correlation in our 
sample of galaxies; a non-Spearman rank order coefficient 
shows a probability of correlation lower than 40 per cent. 

On the other hand, cosmological and semi-cosmological 
chemodynamical simulations by Gibson et al. (2007) and 
Kobayashi (2004) respectively predict slightly shallower val- 
ues, although still steeper than the values obtained here. 
In particular, Gibson et al. (2007) quoted, for a galaxy with 
a ~250 km s _1 a value of grad[E/Fe]=+0.1, while Kobayashi 
(2004) gives mean values of grad[E/Fe]=+0.15 and +0.2 for 
non-major mergers and major mergers galaxies respectively. 

There have been two previous studies presenting [E/Fe] 
gradients with high quality data in two individual galaxies. 
Mendez et al. (2005) found a positive gradient in the galaxy 
NGC 4697 which was reproduced by above quoted models 
of Pipino et al. (2006), and Proctor et al. (2005) analysed 
the gradient of the galaxy NGC 821, which showed a young 
population in its centre, finding a [E/Fe] gradient compati- 
ble with being zero or slightly negative. However, these au- 
thors claimed that, if, as happens in the solar neighborhood, 
oxygen does not track Mg in early-type galaxies, then the 
derived [E/Fe] could turn out to be positive. 

Also Mehlert et al. (2003) derived [E/Fe] gradients for a 
sample of 35 early-type galaxies in the Coma cluster. They 
obtained a mean value of 0.05 ± 0.05 dex, and the devia- 
tion from the mean for all their objects can be explained by 
the errors alone. Therefore, they concluded that early-type 
galaxies show a/Fe gradients consistent with zero. This does 
not exclude, however, the presence of negative a/Fe gradi- 
ents for some of the galaxies of their sample. In any case, 
the quality of this sample of gradients is not as high as the 
one in the other aforementioned studies. These authors also 
analysed the relation between the gradients of metallicity 
and [E/Fe] without finding any significant correlation. 

In order to study whether galaxies with positive and 
negative [E/Fe] gradients are intrinsically different, we have 
analysed the relations between the slope of the gradient and 
other parameters of the galaxies. Fig. [18] shows the relation 

5 Other scenarios can produce differences in the [E/Fe] ratio, as 
differences in the initial mass function or selective loss of gas. 
See, e.g. Faber et al. 1992, Worthey 1998; Trager et al. 2000b for 
details. 



between the [E/Fe] gradients and the central velocity dis- 
persion. A non-parametric rank-order Spearman test gives 
a non-significant correlation, although galaxies with nega- 
tive gradients of [E/Fe] tend to be in the low mass end of 
our sample. We showed earlier in this section that there 
exists a correlation between the [Z/H] gradient and the cen- 
tral value of [Z/H]. In Fig. \T§\ we check the same relation for 
the values of [E/Fe]. With different symbols galaxies with 
power-law (stars) and core (circles) inner profiles are repre- 
sented. This cut also corresponds to a cut in the central a at 
200 km s" 1 . As can be seen, there is a correlation between 
the central value of and the gradients in [E/Fe], but there 
seems to exist a dichotomy between galaxies with core (or 
galaxies with a > 200 km s _1 in our sample) and power-law 
(all galaxies with a < 200 km s _1 in our sample) inner pro- 
files. For galaxies with a power-law inner profile we found 
a clear correlation between the central value and the gra- 
dient of the [E/Fe], while galaxies with a core profile lie 
above the relation. This different relation may be indicating 
fundamental differences in the formation processes of galax- 
ies with a lower and above 200 km s" 1 . This value could 
mark the transition between wet and dry mergers (Faber 
et al. 2005), and its physical motivation may be related to 
the thermal properties of inflowing gas in these galaxies and 
their interplay with feedback processes (Binney 2004; Dekel 
& Birnboim 2006). The transition between the cold flows 
and hot flows can be very sharp, especially if feedback from 
active galaxy nuclei is included. 

When this mechanism is introduced in cosmological 
models of galaxy formation, most massive galaxies form 
their stars at high redshift and in very short time-scales 
and then have assembled later mainly through dry-mergers 
(Cattaneo et al. 2006). This scenario can also explain the 
preferentially boxy isophotes of these systems, their lower 
rotation and their older stellar populations. If confirmed, it 
also supports the idea of biased merging, for which more 
massive galaxies merge, preferentially , with other massive 
galaxies, as predicted in the hierarchical models of galaxy 
formation (e.g. Kauffmann & Chariot 1998) and also sup- 
ported in the relation of the metal poor globular-cluster col- 
ors and the galaxy luminosity (Brodie & Strader 2006). 

We explore now the correlation between the [E/Fe] gra- 
dient and the a4 parameter, which measures the degree of 
boxiness or disciness of the isophote shapes. As can be seen 
in Fig. 1201 the correlation is not as clear as in the case of the 
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Figure 19. Central [E/Fc] inside a r c f/8 aperture as a function of 
the [E/Fe] gradient. Symbols are the same as in Fig. 1101 Stars and 
circles represent galaxies with power-law and core inner profiles 
respectively. 
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Figure 20. Relation between the [E/Fe] gradient and the param- 
eter (a4/a)xl00, which indicates the deviation of the isophotes 
from perfect ellipses. 



[Z/H] gradient. If we exclude the galaxy NGC 2865 there is 
perhaps a trend for which galaxies with a more positive a4 
show, also, a steeper, negative [E/Fe] gradient, but, clearly, 
NGC 2865 does not fit into this trend. This may represent 
a temporal state, as this galaxy shows a very young popula- 
tion in its centre (this galaxy also deviates from the Mgb-er 
relation), but we cannot conclude anything without a larger 
sample. 

Finally, we analyse the relation between the [E/Fe] gra- 
dient and the anisotropy parameter as we did in Sec. 17.21 
Fig. [21] shows this relation. Contrary to that seen with the 
metallicity gradient, there is no significant correlation be- 
tween these two parameters for the whole sample. But if we 
eliminate the two most extreme galaxies (NGC 1600 and 
NGC 2865) a non-parametric rank-order test gives a proba- 
bility of no correlation lower than 1 per cent. 
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Figure 21. Relation between the [E/Fe] gradient and the 
anisotropy parameter for our sample of galaxies. 



Galaxy 


t 




a 


NGC 1600 


5.3 


2.5E 


- 06* 


NGC 1700 


-0.18 




0.429 


NGC 2865 


0.25 




0.400 


NGC 3377 


6.26 


2.9E 


- 09* 


NGC 3379 


21.56 


0.000* 


NGC 3384 


7.21 


2.2E 


- 11* 


NGC 4387 


2.00 




0.025 


NGC 4458 


-0.11 




0.456 


NGC 4464 


4.98 


5.1E 


- 06* 


NGC 4472 


14.84 


1.9E 


- 32* 


NGC 4551 


-2.64 


0.005* 



Table 6. Local metallicity and velocity dispersion correlation, 
t-parameter and significance of the correlation between the local 
[Z/H] and the local velocity dispersion. A value of a < 0.01 in- 
dicates that the probability that the correlation is by chance is 
lower than 1 per cent, which correspond to a 3-cr correlation. 



8 RELATION BETWEEN THE LOCAL 

STELLAR POPULATION PARAMETERS 
AND THE LOCAL VELOCITY DISPERSION 

If galactic winds were the only mechanism responsible for 
the presence of metallicity gradients we might expect to 
find a correlation between the local metallicity and the local 
potential well, and between the metallicity and the [E/Fe] 
gradient. We showed in Sec. [7] that the latter is not present 
in our sample, but several authors have found correlations 
between the colors and metallicity gradients, and the local 
potential gradient, parametrized either using the escape ve- 
locity (Franx & Illingworth 1990, Davies et al. 1993) or the 
velocity dispersion (Mehlert et al. 2003). It is interesting to 
check if this correlation is present in our sample. 

Fig. 1221 shows the relation between [Z/H] at different 
galactocentric distances versus the velocity dispersion mea- 
sured at the same location. To study the degree of corre- 
lation we performed a non-parametric Spearman-rank test. 
Table [6] shows the results of this test. Statistically signifi- 
cant correlations are marked with an asterisk. We have also 
drawn the result of a unweighted linear-fit between the two 
parameters. 

We do not find a strong correlation between the local 
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Figure 22. Relation between the local [Z/H] and the local velocity dispersion. Lines represent the best-lincar-fitting to the data, when 
a significant correlation exists. 



metallicity and the local a for all the galaxies in the sample. 
However, we do find correlations for 7 out of 11 galaxies. 



The correlation between the local metallicity and local 
a arises naturally in a scenario where the star formation pro- 
ceeds until the energy release by the supernova overcomes 
the binding energy and then the gas is expelled from the 
galaxy preventing more star formation. As pointed out by 
Franx & Illingworth (1990), dissipative models with inward 
flows of pre-enriched gas would in principle, tend to destroy 
the correlation. However, Davies et al. (1993) also showed 
that due to anisotropy and rotation effects, the velocity dis- 
persion is a poor indicator of the escape velocity, which could 
be also the reason for the lack of correlation in some of our 
galaxies. 



9 DISCUSSION 

Although almost all galaxies show radial abundance gradi- 
ents, the origin remains a matter of debate. A radial varia- 
tion of star formation rate (SFR), or the existence of radial 
gas flows, or a combination of these processes, can lead to 
abundance gradients in discs (e.g. Lacey & Fall 1985; Koep- 
pen 1994; Edmunds & Greenhow 1995; Tsujimoto et al. 
1995; Chiappini, Matteucci & Gratton 1997). Mergers are 
more complicated, because they depend on several param- 
eters as the mass fractions of the systems merging, the 
amount of dissipation and possible associated star forma- 
tion. However, under different scenarios we would expect 
differences in the strength of the gradients and in their re- 
lationship with other global galaxy parameters. 

We have found, in the present paper, that early-type 
galaxies show null or very shallow age gradients, negative 
metallicity gradients, ranging from grad[Z/H]= —0.19 to 
—0.51 dex, and both, positive and negative, but very shal- 
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low, [E/Fe] gradients (from —0.09 to 0.06 dex). The exis- 
tence of both, positive and negative [E/Fe] gradients rule out 
simple outside-in scenarios where the gradient is an exclu- 
sive consequence of the delayed onset of the galactic winds in 
the central parts as suggested in some studies (e.g., Franx & 
Illingworth 1990; Martinelli et al. 1998; Pipino et al. 2006). 
This is also confirmed by the lack of correlation between 
[Z/H] and [E/Fe] gradients. 

In the present paper, we have found a correlation be- 
tween the metallicity gradient and both the shape of the 
isophotes, and the degree of rotational support. These trends 
are difficult to explain in monolithic scenarios of galaxy for- 
mation, but are well explained in scenarios of mergers where 
the degree of dissipation decrease with the mass. In this 
context, and assuming that the interaction triggers star for- 
mation in the centre of the remnant, the strength and the 
rapidity of the central burst (in general, the degree of dissi- 
pation during the interaction) would determine the shape of 
the isophotes (a4 becomes larger with the degree of dissipa- 
tion, or with less efficient star formation), the central value 
of [E/Fe] (also higher for more efficient star formation) and 
the strength of the [Z/H] and [E/Fe] gradients - the later 
could become negative if the gas is transformed into stars 
very efficiently in very short time-scales (see Thomas et al. 
1999). 

The transition between mergers with gas and com- 
pletely dissipationless mergers is believed to occur at a crit- 
ical stellar mass of 3x10 10 Mq, which is the mass separat- 
ing the red and blue sequence of galaxies (Kauffmann et al. 
2003) in the color-magnitude diagram. This mass corre- 
sponds roughly to the mass of the galaxies showing stronger 
gradients (Vader 1988; Kormendy & Djorgovski 1989; this 
study). The dichotomy observed in the grad[E/Fe]-central 
[E/Fe] plane suggests that the transition is not gradual. In 
order to explain the existence of red and blue sequences 
in the colour-magnitude diagram of galaxies (Kauffmann 
et al. 2003) , cosmological simulations also need to introduce 
a mechanism that sharply terminates the star formation in 
galaxies with stellar masses M> 3x10 Mq. Feedback from 
active galactic nuclei (AGNs) is the most promising can- 
didate for this extra source of heating (Tabor & Binney 
1993; Tucker & David 1997; Granato et al. 2004; Di Matteo, 
Springel & Hernquist 2005; Kawata & Gibson 2005). 

We have also find, in the present study, a correlation be- 
tween the metallicity gradient and both the central values 
of age (for galaxies younger than ~ 10 Gyr) and metallic- 
ity. We have shown that these correlations are not entirely 
due to the correlation of the errors. If we assume, as sug- 
gested above, that abundance gradients are modified by the 
occurrence of major mergers, the fact that these gradients 
are correlated with the central values, and with structural 
parameters as the shape of the isophotes, suggests that the 
young population found in the centre of some of our galaxies 
is the consequence of star formation triggered during major 
mergers. 

This scenario, however, is not free from problems. 
Thomas et al. (1999) and Pipino & Matteucci (2004) have 
shown that is difficult to reach the high values of [Mg/Fe] 
in the centres of elliptical galaxies when massive secondary 
bursts of star formation are superimposed on an old popula- 
tion, if the new stars form from gas with abundances match- 
ing the ones of present day spiral galaxies, unless those stars 



form with an IMF flatter than that of the Salpeter (1955). 
These studies are potentially useful to constraint the amount 
of new stars that can be formed under this scenario. In a fol- 
lowing paper, we will present this information and test the 
viability of this scenario from the chemical point of view. 

The absence of age gradients in most of the galaxies 
also constraints the amount of new stars that can be formed 
in these episodes. However, due to the age-metallicity de- 
generacy, a burst of star formation older than 2 Gyr and 
more metal rich than the underlying population may pro- 
duce very flat gradients in H/3 and stronger differences in 
the metallic indices. For example, the difference between two 
single stellar populations with ages 5 and 12 Gyr, metallic- 
ities [Z/H] — —0.5 and+0.2 and the same [E/Fe] would be, 
-0.024A in H/3, 1.2A in Mgb and 4.lA in C4668 which is 
compatible with the differences between the indices in the 
centre and at one effective radius for the galaxy NGC 3384, 
which shows a null age gradient. 

The amount of gas present during the interaction to pro- 
duce the observed trends and properties of the less massive 
galaxies does not need to be large. Numerical simulations of 
mergers of disc galaxies with a gas fraction for the progenitor 
galaxies of only 10 per cent can reproduce the kinematical 
and photometric properties of intermediate mass ellipticals 
(Jesseit et al. 2006). This is approximately the mass fraction 
of gas remaining in a galaxy at z~0.6 before experiencing its 
last major merger in the cosmological simulations by Meza 
et al. (2003). 



10 SUMMARY 

We present measurements of stellar population parameter 
gradients with unprecedented quality for a sample of 11 
early- type galaxies in the field and the Virgo cluster. 

For most of the galaxies in our sample we find null or 
shallow age gradients, negative [Z/H] gradients ranging from 
(—0.135 and —0.514) and both, negative and positive, al- 
though shallow, [E/Fe] gradients. 

We do not find any strong correlation between the 
metallicity gradient and the central velocity dispersion of 
the galaxies. Galaxies with stronger metallicity gradients are 
those with a ~ 200 km s _1 . For galaxies with a below this 
limit there is a relation for which more massive galaxies show 
steeper metallicity gradients. This trend seems to change at 
a ~ 200 km s _1 , as previously noted by other authors (e.g. 
Vader 1998; Peletier et al. 1990), although we do not have 
enough galaxies with a above this limit to confirm the result. 

We found that metallicity gradient correlates with both 
central metallicity and central age, although the correla- 
tion with the central age is only visible for galaxies with 
mean ages below ~10 Gyr. Galaxies with a younger central 
equivalent-SSP age show also steeper metallicity gradients 
and higher central metallicity. 

Both, [Z/H] and [E/Fe] gradients show strong corre- 
lations with the shape of the isophotes. The shape of the 
isophotes in a merger can be determined by the degree of dis- 
sipation during the interaction (e.g., Bekki & Shioya 1997). 
[Z/H] gradients do not correlate with [E/Fe] gradients, which 
challenges the scenario where the local potential-well is the 
only responsible agent for the presence of gradients in early- 
type galaxies. 
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[E/Fe] gradients show a correlation with the central val- 
ues of [E/Fe], but two different sequences are defined for 
galaxies with central a less and greater than 200 km s _1 . 
Most massive galaxies show a higher central [E/Fe] for a 
given [E/Fe] gradient than less massive ones. 

In general, the presented gradients and their corre- 
lation with other parameters are compatible with a sce- 
nario where elliptical galaxies formed through the merger 
of smaller structures and the relative amount of dissipation 
experienced by the baryonic mass component along ellipti- 
cals stellar mass assembly decreases with increasing galaxy 
mass. However, two results, the trend between the metallic- 
ity gradients and the central a and the relation between the 
gradient and central [E/Fe], suggests that this transition is 
not gradual but there is a sharp change in the properties of 
galaxies above and below a ~ 200 km s _1 . 

Although the quality of the gradients analysed here is 
very high, the sample is not very large. Therefore, the con- 
clusions presented in this paper will need to be confirmed 
with larger samples. 
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APPENDIX A: PROPERTIES OF INDIVIDUAL 
GALAXIES 

NGC 1600: This galaxy is a X-bright field elliptical E3 
galaxy. It does not show any sign of morphological dis- 
turbances (Michard & Prugniel 2004). The shape of the 
isophotes is boxy and has a core shaped inner profile (Ben- 
der et al. 1989; Faber et al. 1997). This galaxy has been 
proposed as a prototype of a merger where the effects of gas 
has not been very important (Matthias & Gerhard 1999). 

NGC 1700: This galaxy shows clear optical signatures 
of a past merger, including shells and prominent boxiness 
(Franx et al. 1989; Whitmore et al. 1997; Brown et al. 2000; 
Lauer et al. 2005). Photometric fine structure at large radii 
(Schweizer & Seitzer 1992) is also indicative of a merger, 
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a velocity reversal ~50" northeast suggest a major event. 
This galaxy shows a counter rotating core in its centre, but 
radially increasing prograde rotating in the main body of 
the galaxy implies that this was not the same event respon- 
sible for the counter-rotating core (Statler, Smecker-Hane 
& Cecil 1996). The strong rotation at large radius and the 
nearly oblate shape are consistent with N-body simulations 
of group mergers (Weil 1995). The disturbances are com- 
patible with a merger of 3 or more stellar systems 2-4 Gyr 
ago. Despite this process, this galaxy follow the same Mg-<r 
relation that the rest of elliptical galaxies (e.g. Trager et al. 
2000a; Sanchez-Blazquez et al. 2006a). NGC 1700 is situated 
in an small group of galaxies with 3 members (Garcia 1993) . 

NGC 2865: This galaxy also shows strong morpholog- 
ical peculiarities (02 = 10.6, Michard 2005), which can be 
indicative of recent interactions. It also shows very young 
mean age in its centre (Hau, Carter & Balcells 1999; Michard 
2005, see Seel}. 

NGC 3377: This galaxy, situated in the Leo I group, 
shows presence of dust (Lauer et al. 2005) and a stellar disc. 
The galaxy has discy isophotes in the inner regions, but 
boxy in the outer regions (Peletier et al. 1990) , and it has a 
power-law central luminosity profile (Faber et al. 1997; Rest 
et al. 2001). 

NGC 3779: This galaxy from the Leo I group, does not 
show any sign of morphological disturbance (Schweizer & 
Seitzer 1992). There is, however, a small nuclear nuclear 
dust ring (van Dokkum & Franx 1995; Lauer et al. 2005), 
and some ionized gas that extends to a radius of 8 arcsec 
(Macchetto et al. 1996). All the globular clusters seems to be 
old and show a broad dispersion in metallicity (Pierce et al. 
2006). The dynamics of this galaxy is explained assuming 
a merger where the gas dissipation has not been important 
(Gebhardt et al. 2000). 

NGC 3384: The galaxy, also from the Leo I group, is 
the only SO of our sample. It shows a circumnuclear disc 
(Sil'chenko et al. 2003) and a disc or a secondary bar aligned 
with the major axis. This galaxy also has a non-symmetric 
distribution of stars in the disc (Busarello et al. 1996). 
NGC 3384 has been proposed as a candidate for a pseu- 
dobulge (Pinkney et al. 2003). 

NGC 4387: This Virgo galaxy shows and old central age 
(~ 12.5 Gyr, Yamada et al. 2006; Sanchez-Blazquez et al. 
2006b). It is a boxy galaxy and it shows a central velocity 
dispersion drop (Halliday et al. 2001; Emsellem et al. 2004). 
This decrease in a in the central parts of galaxies is predicted 
for remnants of an equal-mass merger of two spirals (Bendo 
& Barnes 2000). 

NGC 4458: This low luminosity galaxy, situated in the 
Virgo cluster, also shows a nuclear stellar disc consistent 
with the presence of a central fast-rotating component in 
the stellar velocity field (Morelli et al. 2004) and a hot kine- 
matically decoupled core. The cold central disc of NGC 4458 
does not show any difference in terms of stellar population 
with the main body of the galaxy (Morelli et al. 2004). The 
galaxy shows deviations from the r 1 ' 4 law in their outskirts 
(Michard 1985; Prugniel, Nieto & Simien 1987; Peletier et al. 
1990), which have been explained as a result of tidal inter- 
action with NGC 4461 and NGC 4486. 

NGC 4464: This galaxy, situated in the Virgo cluster, 
shows a power-law inner profile (Faber et al. 1997; Lauer 
et al. 2005). 



NGC 4472: This giant elliptical, situated in the Virgo 
cluster, has X-ray holes or cavities within radii of ~2 kpc 
which may have been produced during a period of nuclear 
activity that began 1.2xl0 7 years ago and may be ongoing 
(Biller et al. 2004). The globular clusters of this galaxy seem 
to be old and coeval but they show a range in metallicity 
from -1.6<[Fe/H]<0 dex (Beasley et al. 2002). 

NGC 4551: This galaxy, one of the members of the Virgo 
cluster, present a power-law inner profile (Faber et al. 1997; 
Lauer et al. 2005). A central decrease in the velocity disper- 
sion has been also found (Halliday et al. 2001). 



APPENDIX B: TRANSFORMATION TO THE 
LICK SYSTEM 

In Fig. IB1I a comparison between the original Lick/IDS in- 
dex measurements of galaxies in common with our data is 
shown. Indices in the Lick/IDS galaxies have been measured 
within an aperture of 4 x 1.5 arcsec (Trager et al. 1998) 
while we tried to match this aperture extracting the spectra 
within the central 4". Therefore, obtaining an aperture of 
4 x 2arcsec. Fig. IB2l show the comparison between the fully 
calibrated data by Sanchez-Blazquez (2004) and this work. 
We selected, from the first study, the measurements within 
an aperture of 1.4 x 4 arcsec. 



APPENDIX C: COMPARISON WITH OTHER 
STUDIES 

Some of the galaxies analysed in this study have been al- 
ready studied by other authors. In particular, we have 5 
galaxies in common with Kuntschner et al. (2006) and 2 
galaxies in common with Fisher, Sadler & Peletier (1996). 
The data by Kuntschner et al. (2006) were acquired with 
an Integral Field Unit (SAURON), but the authors kindly 
provided us with the extracted gradients from their data as- 
suming an slit width and position angle as the one in our 
study. Figures IC1I and IC2I show this comparison. In some 
cases, we have adjusted the zero point as we are interested 
in the gradients. 



APPENDIX D: CENTRAL LINE-STRENGTH 
INDICES 

Table [DTI shows the indices and velocity dispersion measured 
in the spectra extracted within an aperture of 2" x r c ff/8. 
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Figure Bl. Comparison of the line-strength indices for the galaxies in common between Trager et al. (1998) (T98) and this study. Solid 
line shows the 1:1 line while the dashed line shows the calculated mean offset between both samples (see text for details. 
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Figure B2. Comparison of the line-strength indices for the galaxies in common with PSB04 and this study. Solid line show the 1:1 
relation while the dashed line shows the mean offset between samples (see text for details). 
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Figure CI. Comparison of the line-strength indices with radius measured in Kuntschner et al. (2006) (stars) and in this work (small 
open circles). 
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Table Dl. Line-strength indices and velocity dispersion within an aperture of 2" X r e fj/8. 
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0.030 


0.032 


0.051 


0.018 


0.056 


0.010 


0.010 


0.025 


0.031 


0.043 


NGC 1700 


5.568 


1.595 


3.676 


8.906 


1.929 


5.961 


0.121 


0.262 


4.323 


3.151 


2.957 




0.028 


0.016 


0.021 


0.037 


0.013 


0.038 


0.010 


0.010 


0.016 


0.021 


0.023 


NGC 2865 


3.785 


1.267 


3.217 


6.194 


3.189 


5.765 


0.072 


0.177 


3.062 


2.703 


2.679 




0.028 


0.017 


0.022 


0.038 


0.013 


0.040 


0.009 


0.009 


0.017 


0.021 


0.025 


NGC 3377 


5.049 


1.464 


3.413 


7.597 


1.860 


5.531 


0.134 


0.269 


4.501 


2.976 


2.841 




0.018 


0.011 


0.014 


0.029 


0.008 


0.030 


0.010 


0.010 


0.011 


0.016 


0.016 


NGC 3379 


5.522 


1.575 


3.530 


7.961 


1.429 


5.641 


0.157 


0.306 


5.011 


3.057 


2.950 




0.018 


0.012 


0.014 


0.029 


0.008 


0.030 


0.010 


0.010 


0.011 


0.016 


0.016 


NGC 3384 


6.040 


1.709 


3.678 


9.252 


1.910 


6.423 


0.144 


0.287 


4.592 


3.371 


3.307 




0.017 


0.010 


0.013 


0.028 


0.008 


0.029 


0.010 


0.010 


0.011 


0.015 


0.014 


NGC 4387 


5.455 


1.498 


3.361 


6.025 


1.715 


5.480 


0.109 


0.244 


4.224 


3.028 


2.893 




0.032 


0.018 


0.024 


0.043 


0.016 


0.043 


0.010 


0.010 


0.018 


0.023 


0.026 


NGC 4458 


4.855 


1.350 


3.169 


5.154 


1.674 


5.144 


0.100 


0.232 


4.233 


2.730 


2.445 




0.036 


0.021 


0.028 


0.048 


0.017 


0.046 


0.010 


0.010 


0.020 


0.024 


0.028 


NGC 4464 


5.015 


1.460 


3.327 


5.225 


1.595 


5.322 


0.109 


0.248 


4.451 


2.714 


2.612 




0.030 


0.018 


0.024 


0.042 


0.015 


0.042 


0.010 


0.010 


0.017 


0.022 


0.025 


NGC 4472 


5.675 


1.586 


3.557 


8.424 


1.507 


6.081 


0.157 


0.312 


5.064 


3.129 


3.018 




0.020 


0.014 


0.015 


0.030 


0.009 


0.033 


0.011 


0.010 


0.012 


0.018 


0.018 


NGC 4551 


5.691 


1.589 


3.435 


6.845 


1.807 


5.870 


0.114 


0.256 


4.449 


3.150 


3.074 




0.032 


0.018 


0.025 


0.042 


0.016 


0.042 


0.010 


0.010 


0.018 


0.023 


0.025 



Table Dl. Continuation 



